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Ames Research Center 
Moffett Field, Calif. 

SUMMARY 


The response of the magnetosphere boundary in a steady solar wind to 
small initial departures from equilibrium is investigated in accordance with 
the classical model of Chapman and Ferraro. If the wavelength and amplitude 
are sufficiently small that curvature and second -order effects can be disre- 
garded, all perturbations, except those having wave fronts alined with the 
direction of the local magnetic field, are found to damp exponentially with 
time and to drift along the boundary with the tangential component of the 
solar wind. Alined waves, which neither damp nor amplify in this approxima- 
tion, are examined further by inclusion of curvature and higher order effects. 
A first-order analysis shows that curvature introduces a destabilizing effect 
in small regions in the vicinity of the neutral points and a stabilizing 
effect elsewhere. Possible geophysical consequences, such as the persistent 
magnetic agitation of the polar regions, are discussed. An exact solution for 
an alined cylindrical solitary wave having an initial form of a circular arc 
is also presented to illustrate a mode of response that appears to permit 
injection of elongated and widely separated columns of solar wind plasma into 
the magnetosphere under certain conditions and to provide a mechanism for 
momentum transfer from the solar wind to the ambient magnetosphere plasma. 


INTRODUCTION 


This paper is concerned with a theoretical determination of the dynamical 
response of the boundary of the geomagnetic field in a steady solar wind to a 
small initial disturbance from the equilibrium configuration. It differs from 
previous related studies by Dungey (refs. 1 and 2) and Parker (ref. 3 ) in that 
(a) the analysis is based on strict application of the classical theory of 
Chapman and Ferraro and does not invoke additional assumptions, such as those 
brought over from the study of somewhat analogous hydrodynamic problems 
related to the generation of surface waves by wind, and (b) effects of non- 
linear terms and of curvature of the equilibrium shape of the boundary are 
included in those cases that are neutrally stable in the first approximation. 

Although the Chapman -Ferraro theory is itself a highly simplified theory, 
it has emerged from numerous comparisons with experiments conducted in space 
as the leading theory for explaining many features of the interaction between 
the geomagnetic field and the solar wind. It is true that such comparisons 



and other considerations have shown that additional features, such as a 
detached bow wave or viscous or turbulent mixing at the boundary, are probably 
important, but they appear more as additions than alterations to the basic 
theory (e.g., see refs. 4 and 5 for recent resumes). Inasmuch as the present 
application of the Chapman -Ferraro theory is intermediate between the two 
principal previous applications to the interaction .of the geomagnetic field 
with a steady solar wind and with the front of an immense and rapidly advanc- 
ing cloud of solar plasma (see ref, 6 for an extensive review of these appli- 
cations ), it is anticipated that the results will have the same general degree 
of validity and usefulness. At very least, the results have merit in their 
own right as part of a consistent set of formal deductions and consequences 
of the basic Chapman -Ferraro model of the interaction of the geomagnetic field 
and solar plasma. 

The results themselves have a number of features in common with those 
given previously by Dungey and Parker, and some features that differ substan- 
tially therefrom. In particular, the conclusion that the boundary of the 
geomagnetic field is unstable is not supported in general by the present 
results. The results thus appear to be more in keeping with the suggestion of 
Dessler (refs. 7 and 8) that magnetometer and plasma data indicate the bound- 
ary to be stable. (See, however, ref. 9 for a commentary on the earlier of 
these papers.) In some respects, however, the dynamical response determined 
in the present analysis is more complicated than can be described by a simple 
statement of stability or instability. It is found, for instance, that a 
small portion of the boundary over the polar regions is always unstable for a 
certain class of initial disturbances. A major portion of the boundary, how- 
ever, is stable to sufficiently small disturbances, although unstable to 
larger disturbances. The analysis also discloses a mode of response that 
appears to permit injection of elongated and widely spaced columns of solar 
plasma into the magnetosphere under certain conditions and to provide a mech- 
anism for momentum transfer from the solar wind to the plasma contained in the 
magnetosphere . 

A somewhat similar injection mechanism has also been proposed recently by 
Barthel and Sowle (ref. 10 ). Details of the analysis are different in several 
respects, however, and the final conclusion differs in one very Important fea- 
ture. In particular, Barthel and Sowle conclude that injection will occur 
only if there is a very sudden increase in the intensity of the solar wind, 
such as might occur with a frequency of about 10 per year following large 
solar flares. The present analysis suggests, however, that injection will 
occur with even very modest irregularities in the solar wind, such as might 
be capable of forming an elongated dent in the magnetosphere boundary having 
an initial depth of only a few tens of kilometers. 

STATEMENT OF MATHEMATICAL PROBLEM 


The theory of Chapman and Ferraro provides a precise, although idealized, 
representation of the interaction between the geomagnetic field and the solar 
wind, be it steady or unsteady. According to this theory, the solar wind is 
considered to be free from significant influence of magnetic fields, and the 
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earth is considered to he devoid of any atmosphere in the region penetrated 
by the solar plasma. The interaction is such that a current system is estab- 
lished which terminates the geomagnetic field abruptly at a distance of the 
order of ten or more earth radii from the earth, and which leads to forces 
that retard and repel the advance of solar plasma toward the earth. In this 
way, the geomagnetic field is confined to a hollow or cavity carved out of 
the solar plasma. Within this cavity, or magnetosphere, the magnetic field 
satisfies the equations of a static magnetic field in a vacuum, that is 

div B = 0 , curl B = 0 (l) 

The total magnetic field B in the magnetosphere is the sum of the permanent 
magnetic field Bp of the earth and induced magnetic fields B' due to elec- 
tric currents in the ionosphere, the magnetosphere, and at the~ interface 
between the magnetosphere and the solar wind. The permanent magnetic field 
is usually represented in studies of the boundary of the geomagnetic field by 
a simple magnetic dipole placed at the center of the earth and alined with 
the geomagnetic axis as determined from magnetic surveys. The effects of 
electric currents within the ionosphere and the magnetosphere are usually con- 
sidered to be too small, and also too uncertainly known, to be included in the 
analysis, although the effects of the presence of a substantial ring current 
on the location of the boundary of the geomagnetic field in a steady solar 
wind have been investigated by Spreiter and Alksne (refs. 11 and 12). In 
this paper, we are concerned only with the dynamical response of an element 
of the boundary to a small initial displacement from its equilibrium location. 
For this purpose a detailed specification of magnetization and current dis- 
tribution at points far removed from the element is not required. 

A complete mathematical description does require, however, that appro- 
priate conditions be specified at the boundary of the magnetosphere, that is, 
at the magnetopause. Since the location of the magnetopause is unknown and 
must be determined as part of the solution, it is necessary to specify two 
boundary conditions to determine a unique solution. One is that the normal 
component of the geomagnetic field vanish at the magnetopause, that is 

B s . fi s = B n = 0 (2) 

where B s refers to the value of B at the magnetopause and n s refers to 
the unit normal taken positive when directed into the magnetosphere. The 
second is that the pressure exerted on the boundary by the particles of the 
solar wind be balanced by the magnetic pressure B s 2 /8jt. This relation may 
be written as 

B g 2 /8« = Kmn[(V - v) • n s ] 2 (3) 

where m, n, and y refer to the mass, number density, and velocity of the 
ions (generally considered to be principally protons) of the undisturbed 
incident stream, v refers to the local velocity of an element of the magneto- 
pause, and K is a constant coefficient. It is required, in addition that 

(v - v) • n s > 0 (4) 
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in order to avoid the occurrence of regions of the magnetopause that are 
shielded from direct impact by particles of the solar wind, and for which 
equation (3) would be inappropriate. 

It has historically been considered that the particles are, in effect, 
specularly reflected at the boundary in which case K = 2 (see refs. 13 and 14 
for a discussion of other values for K that have sometimes been used for this 
case). It has been argued recently by Axford (ref. 15), Kellogg (ref. l6), 
Spreiter and Jones (ref. 17), and others, however, that the presence of a weak 
interplanetary field may have the important effect of causing the solar wind 
to behave as a continuous fluid over length scales that are large compared 
with the proton Larmor radius. Since the Larmor radius is about 600 km for a 
300 km/sec proton in a magnetic field of 5 7 and a typical dimension of the 
magnetosphere is of the order of 10 5 km, it follows that the solar wind would 
flow around the magnetosphere much as a fluid about a solid object. Since, 
moreover, the effective Mach number based on the speed of a magnetoacoustic 
wave is much greater than unity, there would be a detached bow wave upstream 
from the magnetosphere. Nevertheless, the aerodynamic pressure exerted by the 
wind on the boundary would still be represented, according to the Newtonian 
theory of hypersonic flow, by the right side of equation (3) with K = 1. 

This change from K = 2 to K = 1, without further change of the governing 
equations (l) through (4), leaves the calculated form of the magnetopause 
unaltered, although all linear dimensions from the center of the earth to the 
magnetopause increased by a factor of 2 1//s . Results so obtained have been 
shown (IG Bulletin 84, June 1964) to be in excellent agreement with data 
obtained in space with IMP -I satellite (Explorer XVIII). The effects of such 
a change are of even less importance in the present investigation, since we 
are concerned with deviations from the equilibrium configuration, and K 
cancels out of the analysis in most cases. 


The equilibrium steady-state configuration of the magnetopause and the 
enclosed magnetic field is defined by the same set of equations, except that 
v is equated to zero in equations (3) and (4). Equation (3) is thus replaced 




B| e /8jt - Kmn( V . n se ) 2 


(5) 


y 



Figure 1.- View of coordinate system and an 
element of the magnetopause. 


where the subscript e indicates 
values associated with the equilibrium 
configuration. 

Solutions of the foregoing equa- 
tions to be presented in the remainder 
of this paper are concerned with the 
behavior of limited portion of the mag- 
netopause, and are expressed in terms 
of the local rectangular Cartesian 
coordinate system illustrated in fig- 
ure 1. In all cases, the origin is 
placed at the center of the region of 
interest, either at or near the mag- 
netopause, and the z axis is directed 
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normal to the boundary and into the magnetosphere. It is thus convenient to 
represent the coordinates of the boundary by z = Z(x,y,t) where t represents 
time. Further, attention is confined to a sufficiently small portion of the 
magnetopause that Z is small compared not only with the over -all dimensions 
of the magnetosphere, such as the distance from the center of the earth to the 
nearest part of the magnetopause, say 60,000 km under representative condi- 
tions, but also with the dimensions in the x and y directions of the mag- 
netopause element under consideration. The velocity of the undisturbed 
incident solar plasma is expressed in terms of the components V x , V y , V z , and 
the unit vectors 2, 3> k parallel the x, y, z axes, respectively; thus, 

V = iV x + jV y + kV z (6) 

The equilibrium magnetic field can be expressed similarly as 

He = iBe x + 3B ey + kBe z (7) 

where Be x = B e cos 0, ®e y = ®e sin ® > and Be z » 0. The Cartesian components 
of the unit normal ft s to the magnetopause are 

- _ i(-dZ/dx) + 3(-dZ/dy) + k 
[l + (dz/dx) 2 + (dZ/dy) 2 ] 1/ ' 2 

Since, furthermore, 

v • h s = (dz/dt)k • h s 


( 8 ) 


(9) 


at any point xy on the magnetopause, 


-V^Z/dx - V az /by + v z - bZ/bt 

(V - v) ' n q = * 7 — (10) 

[1 + (bZ/bxf + (bZ/by) 2 } 


The foregoing equations completely describe the mathematical model on which 
the remainder of -the discussion is based. The remainder of this paper is 
concerned with the determination and discussion of solutions for the dynami- 
cal response of the magnetosphere boundary to a variety of initial displace- 
ments from equilibrium. It is assumed throughout that the undisturbed 
incident solar wind is uniform and steady in space and time. 


FIRST-ORDER SOLUTION FOR SMALL AMPLITUDE WAVES ON AN 
OTHERWISE FLAT MAGNETOSPHERE BOUNDARY 


Consider first the case of small amplitude waves of sufficiently short 
wavelength that curvature of the equilibrium magnetosphere boundary can be 
disregarded. Neglecting the effects of curvature is permissible under most 
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circumstances in which the wavelength is small compared with the radius of 
curvature of the equilibrium boundary. Under representative conditions, the 
latter dimension is of the order of 60,000 km or more everywhere except in a 
small region in the vicinity of each neutral point. It is further assumed 
that the ratio of the amplitude to the wavelength is sufficiently small that 
we can safely disregard all terms of higher order than the first in Z and 
its derivatives, where z = 0 is taken to represent the coordinates of the 
flat equilibrium boundary, and z = Z(x,y,t) represents those of the perturbed 
boundary. 

Under the above conditions, we can approximate (V - v) • n s by the 
numerator of the right side of equation (10), and its square by 


[(V - v) 





'Yx az az _i_ az\ 

v z Sx v z a y v z at/ 


( 11 ) 


The equilibrium magnetic field associated with the undisturbed flat boundary 
is simply B e = const. The total field can be expressed by 


B = Bg + b 


( 12 ) 


where b is a function of space and time that represents the perturbation 
magnetic field due to the departure of the magnetosphere boundary from its 
assumed flat equilibrium shape. To first order in b, we have 


B 2 = 




+ b 


z 


2 


* ®e x + y + 2^Be x b x + B e yby^) 

= B e 2 [l + 2(b x cos 9 + by sin 9 )/ B e ] (13) 


Since the difference between the value for b at the actual location of the 
boundary and at the equilibrium position (z = 0) is given by an expression of 
the form (db/dz)Z and is, hence, a product of small quantities, it is suffi- 
cient to evaluate B s 2 at z = 0; thus, 

B s 2 = B e 2 [ 1 + 2(b x cos 9 + by sin S ) /B e ] z _ c (14) 
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Substitution of equations (ll) and (l4) into equation (3) yields 


B e 2 

8jt 


1 + 2 


\>x cos 0 


Be 


+ by sin 9 


z=o- 



Since the equilibrium surface lies in the xy plane, n se = k for this case, 
and equation (5) reduces to 


B e 2 /8it = KmnV z 2 


( 16 ) 


Combining equations ( 15 ) and (l6) yields 



cos 0 + by 

B e 


sin 


0 s 

J 


' z=o 


'Vx az Yi. _i_ az\ 

V z 8x + V z 8y + Vz Sty 


(17) 


The right side of equation (17) can be simplified by transforming to a coor- 
dinate system fixed with respect to the tangential component of the solar 
wind. We thus introduce the new variables 


| = x - V x t , n = y - V y t , 


i = z , t = t 


( 18 ) 


Equation ( 17 ) thus becomes 



cos 0 + b^ 
Be 


sin 0 




_i_ 8z 
v z 8t 


( 19 ) 


In order to proceed further, it is necessary to determine the appropriate 
relation between b and Z. Substitution of equation (12) into equation (l), 
with B e = const, shows that div b = 0 and curl b = 0. This implies that 


b = -grad ft 


( 20 ) 


where Q. is a solution of Laplace’s equation 
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v 2 n = o 


(21) 


The boundary conditions for b are that 


(5) 


^ = 00 


= 0 


(22) 


and 



dZ = )Zd[ 
ds dj; ds + 


dZ dri dZ 
br\ ds ” 5f* 


Sz 

cos 0 + >r— sin Q 

OT] 


(23) 


The corresponding conditions for Q are 

(&)<, = const (24) 

f — 00 


which can he equated to zero without loss of generality, and 


JL 

' B e W 5=0 


dz 

= ST 003 e 



sin 0 


(25) 


A general solution of equation (21) that vanishes at infinite £ can he 
found hy separation of variables to he a summation of terms of the form 


2 = (A sin ki| + B cos ki|)(C sin k^ + D cos k2q)exp[ -(kj. 2 + k2 2 ) £] 

( 26 ) 


in which the integration constants may vary with time. If attention is con- 
fined to wavelike boundary perturbations having nodal lines at t = 0 
arranged in the form of a rectangular grid, and the coordinate system is ori- 
ented without further restriction so that the £ and r\ axes are parallel the 
nodal lines at that time, we can write 

(z ) T=0 = e(o) cos ki| cos k^r) ( 27 ) 

The proper evaluation of the coefficients of equation (26) so as to satisfy 
equations (25) and (27) yields the following expression for 0 . at r = 0: 
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(28) 


(fl) 


T=0 


= -e(0) ( ~ cos 


0 sin ki| cos 


k 2 

k 2 r\ + — sin 


0 cos ki£ sin k^n ! e 




where 


k = (k, 2 



(29) 


The left side of equation (19) thus becomes, at t = 0 


f bp cos 0 + b n sin 0\ e (0) r , n 2 2 „ , 2 . 2 ^ „ , 

LJ = -“-[(ki cos^ Q + k 2 sin 0)cos k^ cos k 2 rj 

V B e /^ =0 k 

T=0 

- 2k x k 2 sin 0 cos 0 sin kii sin k 2 r]] (3O) 


It is convenient to restrict attention at this point to the not inconsid- 
erable class of cases in which any one of the following additional conditions 
hold 

sin 0=0, cos 0=0, ki=0, k 2 =0 (31) 

The first two conditions permit the wavelength to be arbitrary in both direc- 
tions but require that the equilibrium magnetic field be parallel to either 
of the two sets of nodal lines of the initial boundary perturbations. The 
latter two conditions permit arbitrary orientation of the nodal lines and the 
equilibrium magnetic field, but require that the boundary perturbation be 
two-dimensional, that is, have the form of cylindrical corrugations. Equa- 
tions (19) and (3O) indicate that the nodal lines remain fixed for all time 
for these cases, and that the solution is given by 


Z = s(t) cos kii cos k 2 r) 


(32) 


where the amplitude e(r) satisfies the equation 


d€ -V 2 (ki g cos 2 0 + k 2 2 sin 2 0)e 

dT k 


-Kq€ 


(33) 


The solution is 


= e n e 


-KqT 


(3^ 


e 
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where e 0 represents the amplitude at t = 0; that is, e Q = e(0). It should 

be observed that the behavior of the boundary is considerably more complicated 

if none of the four conditions given in equations (31) apply, since then 

(SZ/St) fails to vanish on the nodal lines where (Z) = 0 and the form of 

t=o t=o 

the boundary perturbation depart immediately from the original cos k x £ cos k 2 r) 

form. 


Substitution into equations (32), (33)> and (34) of the original physi- 
cal variables using the relations given by equation (19) yields 

-t/t* 

Z = e 0 e 0 [ cos ki(x - V x t) ] [cos k 2 (y - V y t)] (35) 

where the time constant t* is given by 


and 



k 

V z (ki 2 cos 2 0 + k 2 2 sin 2 0) 


(A x 2 + A y 2 ) 1/2 
2*V Z 



A x A y 


. 2 _ 2 
sin 0 +- Ay 



( 36 ) 


A x = 2it/ki , Ay = 2jt/k 2 


(37) 


represent the wavelengths measured in the x and y directions, respectively. 
Two general conclusions follow immediately from the above results. The first 
is that the presence of the quantities x - V x t and y - Vyt in the arguments 
of the cosines of equation (35) indicates that the wave system drifts along 
the magnetosphere boundary with the tangential component of the solar wind. 
The second is that the time constant t* is positive, indicating that the 
waves damp with time, or at least never grow. 

If ki or k 2 vanishes, the boundary perturbations are two-dimensional 
corrugations with the nodal lines and other generators parallel to the x or 
y axis, respectively. The time constants for these cases are given by 


(tS) 


A x =°° 


A, 


2kV z sin 2 0 


A 


x 


(to), - 

Ay— 0° ^ tt 2 


(38) 


2rtV z cos 0 
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A plot of the variation of (toK 

Ay =oc 4 0=75° 

with wavelength Ax for various orien- / 60 * 

tations 0 of the equilibrium magnetic « / / 

field is shown in figure 2 for a repre- *. 3 _ / / 

sentative case in which V z = 500 km/sec. t / ^ 45 ° 

We see, for instance, that the time con- / / s' 

stant for decay of corrugations having 2 . / / s' sbo° 

east-west generators (Ay = oo) and sit- § \ / s' 

uated near the equatorial plane (0 = 0°) | / z' s' ^ °° 

is of the order of a second for wave- “ i - / s^ 

lengths as great as 3000 km, and propor- i v z =soo km/sec 

tionally less for shorter wavelengths. 

It increases to infinity, implying no , 1 — 1 1 1 

damping, as 0 goes to 9 0 u and the ° l0 °° 4 °°° 5 °°° 

nodal lines are rotated so as to become 

parallel to the geomagnetic field lines. T °^ time . con f ant vith 

. 0 0 wavelength for two-dimensional corruga— 

This example illustrates the general tions on a flat equilibrium boundary. 

result that corrugations having genera- 
tors alined with the equilibrium field are neutrally stable to first order in 
small perturbations. This case will, therefore, be re-examined to higher 
order of accuracy subsequently herein. 


If, on the other hand, neither k x nor k 2 vanishes, but 0 = 0° or 90°, 
so that the nodal lines form a rectangular grid and the equilibrium field 
lines are alined with one set of nodal lines, we have the following expres- 
sions for the time constant 




+ VV 


(t °Ws 


* y> 1/e aa 


/Ay/xx = co 


This result shows the time constant is inversely proportional to 

proportional to the mean scale of the wave system (A x 2 + Ay 2 ; 2 
in addition, on the ratio of wavelengths A 

in the x and y directions. A plot of 4 |~ X * / 3 
the results for (t£to , for o J / 2 

o/ 0=jt/2 8 / / 

V z = 500 km/sec is shown in figure 3* $J3- / / / 

We see, once again, that the time con- ^ / / / 

stant is of the order of a second for a 1° / / / 

wide variety of conditions when mean k" 2 ~ / / / 

wavelengths of the order of a few thou- | / , / s' 

sand kilometers are considered. It I / // s' 

should be noted that figure 3 is not Z ' ~ / 

appropriate for displaying the time con- f 

stants for a two-dimensional corrugation ^ 

since then (A x 2 + Ay 2 ) 17 2 goes to inf in- 0 2000 4000 60 2 00 

ity and an indeterminate form appears. mean wavelength, (x x +x 

Results for this case are given correctly Fi& " re 3 ’ ’ v f ia ^°? ° f tl f c 
by equation {39) > however, or, more on a flat equilibrium bounc 


2000 4000 6000 8000 

MEAN WAVELENGTH, (\2+x2) l/2 , km 


Figure 3. - Variation of time constant with 
mean wavelength for rectangular wave grid 
on a flat equilibrium boundary. 
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explicitly and more generally with respect to the relative orientation of the 
directions of the nodal lines and the equilibrium field, hy equations (3&)« 


From the above, we conclude that long wavelength waves die out slower 
than short wavelength waves, that rectangular grid distortions die out slower 
when the wavelength is longer in the direction of the field than in the direc- 
tion normal to the field, and that two-dimensional distortions having nodal 
lines alined with the equilibrium field direction do not damp at all, when 
only first-order terms are retained in the analysis. It follows in addition, 
that, if the boundary is initially disturbed in such a way that it can be 
represented by an arbitrary distribution of waves having the forms described 
above, all components would damp with the passage of time, except those that 
are two-dimensional and have their nodal lines alined with the equilibrium 
field direction. These alone would remain, and they would drift along the 
boundary with a speed equal to the tangential component of the solar wind, 
neither damping nor amplifying as they travel. Since this component vanishes 
at the nose of the magnetosphere near the subsolar point and increases to the 

, o 2 o 1/2 

full speed V -= (V x + Vy + V z ) of the undisturbed incident stream far 
downstream of the earth, the waves tend to drift very slowly near the magne- 
tosphere nose, but accelerate to several hundred kilometers per second along 
the flanks. The results for all the waves considered show, moreover, that 
the time constant is inversely proportional to the normal component V z of 
the solar wind velocity. Since V z equals V at the magnetosphere nose, 
and decreases to zero along the flanks, the time constant grows indefinitely 
with distance along the magnetosphere boundary from the nose, indicating that 
boundary perturbations damp much more slowly on the flanks of the magneto- 



sphere than on che nose. We thus find 
that all except alined waves tend to 
damp rapidly with very little drift 
near the nose of the magnetosphere, 
but damp slowly as they run rapidly 
along the flanks . 

Figure 4 has been included to give 
a more quantitative understanding of 
the degree to which the time constant 
and drift velocity change along the 
equatorial portion of the magnetosphere 
boundary. The solid curve shows the 
equatorial trace of the magnetosphere 
boundary calculated under the usual 
assumption that the geomagnetic field 
in the equatorial plane varies as the 
inverse cube of the distance from the 
center of the earth (Ferraro, ref. l8; 
Beard, ref. 19; Spreiter and Briggs, 


Figure 4.- Variation of lateral, coordinate of 
magi let.o sphere boundary, (r/r,,)cin <y , time 
constant amplification (V/V z ) , and drift 
velocity (V x /V) with geocentric distance 
along a line parallel the solar wind direc- 
tion. 


ref. 14; etc.), and normalized by 
dividing by the geocentric distance 
r 0 to the magnetosphere nose. The 
dashed curves show the time constant 
amplification factor V/V z and the 
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normalized drift velocity V x /V as a function of (r/r Q )cos cp where r rep- 
resents the geocentric distance to the boundary and cp is a longitude angle 
measured from the direction of approach of the undisturbed incident solar 
wind. They are calculated using the relations 


V/V z = (r/r 0 ) 3 

(40) 

« !/ 2 


V x /V = [1 - (r 0 /r) 6 ] 

(41) 


FIRST-ORDER SOLUTION FOR SMALL AMPLITUDE WAVES ON A SLIGHTLY 
CURVED MAGNETOSPHERE BOUNDARY 


The preceding analysis will now be extended to include the effect of 
slight curvature of the equilibrium magnetosphere boundary. We again neglect 
all terms of higher order than the first in Z, where Z(x,y,t) = Z - Z e rep- 
resents the displacement Z of the perturbed magnetopause relative to the 
displacement Z e of the curved equilibrium surface and confine attention to 
the case in which the radius of curvature of the equilibrium surface is much 
greater than both the wavelength and amplitude of the boundary perturbations. 
Under these circumstances, it is permissible to represent the effects of 
boundary displacement from z = 0 at any point by superposition of effects of 
the equilibrium curvature displacement and the perturbation displacement. The 
equilibrium magnetic field can no longer be regarded as a constant in a region 
extending several wavelengths in all directions inside the magnetosphere from 
the origin, as indicated previously in connection with equation (12), however, 
since the irrotationality condition given by equation ( 1 indicates a non- 
vanishing field gradient at a curved boundary. The consequences of this con- 
dition can be illustrated in terms of a local coordinate system in which the 
a axis is alined with the equilibrium magnetic field vector and the v axis 
is the normal to the a axis in the plane determined by the local curvature 
of the field line . Then curl B = 0 evaluated at the equilibrium surface 
leads to 



(te) 


where B S e refers to the magnitude of Be at the equilibrium surface, and 
R refers to the local radius of curvature of the field line. Since v is 
thus very nearly parallel to the direction of z for the conditions described, 
we have the following expression for the equilibrium magnetic field B e in 
the vicinity of the magnetopause 


B 


e = §se 




§se 



( 43 ) 
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The total field at the boundary is thus 


gs 


Sse 




+ be 


(44) 


where b s is the perturbation field associated with the departures Z from 
the equilibrium surface. The perturbation field may be evaluated in the same 
way as in the preceding section for the flat equilibrium boundary. The result 
for a perturbation shape defined by Z(x,y) is found to be the same as deter- 
mined before for a perturbation shape defined by Z(x,y) . The counterpart of 
equation (15) is thus 


B' 


se 

& 


i+ f +2 


b x cos 0 + by sin 9 
B S e 


Z=0 J 


= KmnV^ 


1-2 


Zx az Yz az i_ az 
v z ax v z a y v z at 


(45) 


Subtracting the corresponding expression for the equilibrium state 


B' 


se V TT 2 

ar = KmnV z 


y x az e v v az e 

1 _ 2 ( — - + JL £ 

' v z ax v z ay 


(46) 


and retaining only terms of first order in b x , by, and Z yields 


Z ^bx cos 0 + by sin 9 ' 
5 + V B se 


z=o 


x az yy az i az 


y z ax v z a y v z at 


(47) 


where B se is a constant equal to the value of B se at, say, the origin. 

Introducing the variables |, q, £ , t, defined by equations (l8) trans- 
forms equation (47) into 


Z 

R 


+ 


b^ cos 9 + b^ sin 9 


B, 


se 


Z=0 


_i az 

Vz 8 t 


(48) 


If attention is restricted to the class of cases for which any one of the four 
conditions given by equations ( 31 ) holds, the same line of reasoning as that 
described previously to solve equation ( 19 ) shows that the solution is given 
By 
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Z = € ( t) COS kil COS k2ll 


(^ 9 ) 


where £(t) is given by 


e 


e 0 e 


-T/t* 


(50) 


e 0 represents the amplitude at t = 0 , and 


t* 




k x 2 cos 2 9 + k 2 2 sin 2 

_ 





2jt (Ay 2 cos 2 9 + A x 2 sin 2 6 ) 
A x A y (A x 2 + A /) l/2 




1 + V z t£/R 

(51) 


where A x and Ay represent the wavelengths in the x and y directions as 
defined by equation ( 37 )> and t£ is the time constant defined by equa- 
tion ( 36 ). In terms of the original variables, equations (49) and ( 50 ) are 


— -t /t* 

Z = e Q e ' [cos k x (x - V x t)][cos k 2 (y - Vyt)] 


(52) 


This result shows, again, that the wave system drifts along the magneto- 
sphere boundary with the tangential component of the solar wind, and that the 
amplitude damps with a time constant t*. Now, for R = 00 , t* is always 
positive, and the wave amplitude can never grow, as discussed in the preced- 
ing section. When R is finite, however, more possibilities arise. If R 
is positive (the center of curvature of the equilibrium surface is on the 
magnetic field side of the boundary, that is, inside the magnetosphere), 
t* < t* and the. effect of curvature is to reduce the time constant or, equiv- 
alently, to increase the rate at which the disturbed boundary returns toward 
the equilibrium shape. If, on the other hand, R is negative, the possibil- 
ity exists that t* may become negative, indicating that boundary perturba- 
tions grow in amplitude with time. The critical radius R cr for neutral 
stability is obtained by equating t* to infinity and solving. The result 
is 


R, 


•cr 


= -V z t* = 


A x A y (A x g 


+ A/) 


1/2 


[_2:rt(A x 2 sin 2 9 


+ Ay 2 


cos 2 0)J 


(53) 
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If the nodal lines of the perturbation wave system form a square grid, 
that is, A x = Ay = A, equation (53) reduces to 


R cr = ->/i\/2rt 


(54) 


The formal analysis thus indicates that the boundary is unstable to small 
disturbances if the center of curvature is outside the magnetosphere and the 
radius of curvature of the equilibrium surface is less than \[2k times the 
perturbation wavelength. Since the latter is required to be much less than 
the radius of curvature at the outset of this analysis, it must be concluded 
that all square grid boundary perturbations for which the present results are 
applicable are stable in spite of possible destabilizing influences due to 
curvature of the equilibrium surface. 

If k x or k 2 vanishes and the boundary perturbations are two-dimensional 
corrugations as noted previously, the time constants are given by 


(t*) 


\ x =00 


„ (l 2jt sin 2 d 

L 2 4 v . 


(55) 


(t*) 


Ay=0O 



/ 1 2n cos 2 6 \ 

V R + 4 X j_ 


(56) 


Since A x and Ay are supposed to be very much smaller than |r| in the pres- 
ent analysis, we see that instability (negative t*) occurs only for very 


small 0 if Ax = or for 0 very 


xy 

R 



Figure 5. - Stability boundary for two - 
dimensional corrugations (A x - 00 ) on a 
curved equilibrium surface. 


nearly equal to xc/2 if Ay = oo. For 
example, if A x = oo, A = 1000 km, and 
R = -60,000 km, | sinSj = (l20n)“ 1/2 
= 0 . 0515 - The critical values for 6 
are thus about ±3° for two-dimensional 
boundary perturbations having nodal 
lines alined with the x axis. We 
can conclude alternatively by a simple 
rotation of the coordinate system that 
only those waves having their nodal 
lines within 3° of the direction of 
the field lines will grow with time. 
All others will damp out, and even- 
tually disappear. For other values 
for Ay and R, the numerical results 
will, of course, differ. Figure 5 is 
a plot of the variation with 9 of 
the critical value of the ratio Ay/R 
for waves having A x = 00. Two- 
dimensional corrugations having Ay/R 
less than that indicated by the curve 
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are unstable and, conversely, as indicated. Thus, on a flat surface for which 
R = oo and Ay/ R = 0, all two-dimensional waves are stable, except those per- 
fectly alined with the field direction (0=0°) which have neutral stability. 
All waves on equilibrium surfaces with R > 0 are also stable, but those suf- 
ficiently closely alined with the field are unstable on equilibrium surfaces 
with R < 0. Although R > 0 nearly everywhere on the magnetosphere boundary, 
both cases occur. In particular, negative values for R occur in limited 


regions in the immediate vicinity of the 
figure 6 by the dotted portions of 
traces of the magnetosphere boundary 
(ref. 19)- Also included and indicated 
by dashed lines are several magnetic 
field lines for the noon -midnight merid- 
ian plane, as given by Mead (ref. 20). 

Since only the nearly alined two- 
dimensional waves are significantly 
affected by curvature of the equilibrium 
surface, it is of interest to examine 
the special form to which equations (55) 
and (56) simplify for perfectly alined 
waves, namely 


neutral points N, as illustrated in 



Figure 6. - Traces of the magnetosphere bound- 
ary in the equatorial plane , noon -midnight 
meridian plane, and several intermediate 
planes (ref. 19 )- Magnetic field lines in 
the noon -midnight meridian plane (ref. 20) . 


t* = R/V z (57) 

The time constant is thus equal to the time required for a particle with 
velocity V 2 to travel a distance R. Representative values for conditions 
at the nose of the magnetosphere may be R = 60,000 km and V z - 500 km/sec, 
for which t* = 120 sec. The time constant increases beyond all bounds at 
the neutral points and also far downstream along the flanks of the magneto- 
sphere boundary, since V 2 approaches zero and R does not vanish. Alined 
two-dimensional boundary perturbations thus neither damp nor grow as they 
drift along the magnetosphere boundary in these regions. A short distance 
downstream from the neutral point, we may have R = -60, 000 km and 
Vz = 100 km/sec, in which case t* = -60 sec, and the waves grow with time. 

In all cases, the time constant is long compared with the values representa- 
tive of less elongated rectangular grid perturbations. 

The small regions of instability near the neutral points are of particu- 
lar interest because they provide a means by which solar plasma can be con- 
stantly admitted into the magnetosphere, and may very well be associated with 
the constant magnetic agitation of the polar regions, as summarized, for 
instance, by Fukushima (ref. 21). He assembled a considerable body of ground- 
based magnetometer data gathered during the IGY (July 1957 - December 1958) 
and the Second Polar Year (August 1932 - August 1933 )> and analyzed it in var- 
ious ways according to the time of day and year, and the degree of disturbance 
of the world -wide geomagnetic field. The results show that the polar regions 
are constantly disturbed, even at periods when the planetary geomagnetic field 
is very quiet. Two of the many summary plots presented by Fukushima are 
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reproduced in figure 7- These show the 
latitudinal distribution of the average 
geomagnetic disturbance near the noon 
and midnight meridians for the case of 
Kp = 0 o , that is, for the very quietest 
times as indicated by the planetary 
geomagnetic index, at two selected sets 
of times extending over the entire 
period of the IGY. The results show 
that the magnitude of the polar -cap 
disturbance is of the order of a few 
or several tens of gammas, that both 


Figure 7-- Polar magnetic disturbance near 
the noon and midnight meridians at times 
during the IGY when Kp = 0 o and the 
sun -earth line is inclined at 0° and 
34. 5° to the geomagnetic equatorial 
plane (ref. 21). 


polar caps are equally disturbed when 
the sun -earth line is perpendicular to 
the geomagnetic dipole axis, and that 
the disturbances are a maximum at sum- 
mer noon, and a minimum at winter 
midnight. These results, when taken in conjunction with the corresponding 
changes in the locations of the neutral points, such as calculated by Sp re iter 
and Briggs (ref. 14), appear consistent with what might be anticipated on the 
basis of direct invasion of charged particles or propagation of hydromagnet ic 
waves into the polar regions along the field lines from the vicinity of the 
neutral points on the magnetosphere boundary. 


EXACT SOLUTION FOR ALINED WAVES ON A FLAT 
MAGNETOSPHERE BOUNDARY 


The anlysis presented in the preceding sections is restricted to consid- 
eration of only first-order terms in the perturbation quantities. It is thus 
sufficient for the analysis of decaying waves of small amplitude, and of the 
initial behavior of amplifying waves. It cannot be used, however, to provide 
information on the subsequent behavior of amplifying waves. Although the 
general treatment of such cases presents a very complex problem, the analysis 
of the special case of alined two-dimensional waves on a flat equilibrium 
surface, while nonlinear, remains tractable. This case is, moreover, an 
informative one since such waves were found in the first-order analysis to be 
neutrally stable, and a more accurate analysis is required to elucidate their 
properties . 

For sake of definiteness, consider the equilibrium magnetic field to be 
in the direction of the positive y axis, and that the distortions of the 
boundary surface change with x, but not with y. For such two-dimensional, 
or cylindrical, alined boundary distortions, the perturbation magnetic field 
vanishes and the total field remains constant and equal to the value B e 
for the flat equilibrium surface. Equation ( 3 ) thus becomes 


B e 2 

— n— = Kmn 
Ort 


( -V x dZ/8x + V z - 8Z/8t) 2 

1 + (8z/8x) 2 


(58) 
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and when combined with the corresponding equilibrium relation given by equa- 
tion (l 6) yields 


/. v x 2 \ 

/ dzy i 

/6Z\ 2V X 

6z 

2V X 6z 6z 2 6z 


W/ v z 2 

\av + v z 

Sx 

2. Sx V Z St 

V Z 


Introducing the new variables 

g = X - V x t y T = V Z t , Z = Z - T 


(59) 


(60) 


yields 



(61) 


It can be demonstrated by differentiation that a solution of equation (6l) is 

(Z - A) 2 + (| - B) 2 = (t - C) 2 (62) 

The equivalent expression in terms of the original variables is 

(Z - v z t - A) 2 + (x - V x t - B) 2 = (V z t - c ) 2 ( 63 ) 

This solution represents a circle with center at z c = V z t + A, x c = V x t + B, 
and radius |V z t - C| . We can set B = 0 without loss of generality by locat- 
ing the coordinate system so that the center of the circle is on the 
| = x - V x t = 0 axis and, hence, drifts with the tangential component of the 
solar wind. The remaining integration constants A and C can be evaluated 
by letting t = T at the instant when the center of the circle is at the 
origin, and noting that z = = ±e 0 at | = 0 for all time. The latter 

condition follows from the fact that Sz/6 £ = 0 at 1 = 0 for the choice of 
coordinate system described above, together with the consequence indicated 
by equation (59) that 6z/6t = 0 at such a point. In this way equation (63) 
becomes 

[Z - V z (t - T)] 2 + (x - V x t) 2 = [V z (t - T) - Zj 2 (64) 

The boundary perturbation is thus described by a circular arc of radius R 
given by 
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R = |V z (t - T) - Zm| 


( 65 ) 


with its center located at 


X C = v xt , z c = Vz(t - T) 


( 66 ) 


The center thus moves with the velocity components of the undisturbed incident 
solar wind V x and V z . Although this solution is not useful for determining 
the dynamical response of the magnetosphere boundary to an initial disturbance 
in the form of a periodic wave, it is appropriate for the study of the response 
of the boundary to an initial displacement in the form of a solitary segment 
of a circle. This is not too restricting to render the solution uninterest- 
ing, however, since the region of maximum displacement of any continuous dis- 
tortion will approximate a circle in form. 


Two cases are possible, depending on whether the circular segment bulges 
in the direction of the plasma or the field side of the boundary, or equiva- 
lently on whether Zm is negative or positive. The case of negative 

illustrated in figure 8 will be con- 
sidered first. An expression for the 
half width W at t > T for this case 



can be obtained readily from equa- 
tion ( 63 ) by equating x to W, Z to 0, 
and solving for W. It is 

2 "ix/2 

W - v z (t - T){[l f Vz( t°_ T ) j - l } 1 2 

( 67 ) 

At great time, that is when 
£o/Vz(t - T) « 1, W grows with time 
approximately as 

w * [2V z e 0 (t - T)] l/2 (68) 

With V z = 500 km/sec = 5x10 ^ cm/sec, 

G o = 100 km = 10 y cm, and an initial 
half width Wq = 1000 km = 10 s cm which 
according to equation (68) occurs at 
tq - T = 10 sec, W in cm is 


20 


(69) 


w = 


t - t \ 

ti - i ) 


1 / 2 


XlO £ 



Xio 8 


cm 


It thus requires 30 sec to double the width of the' bulge, and 150 sec to 
quadruple the width. Since the maximum displacement \^a\ remains constant 
as the width increases, the bulge clearly smooths out with the passage of 
time. It should be observed, however, that this process is slow compared with 
the rates charcteristic of unalined waves. 


We next consider the alternative 
case of a dent into the magnetic field 
side of the boundary. Such a boundary 
perturbation is characterized by a pos- 
itive value for = e Q a nd negative 
t - T, as illustrated in figure 9* The 
half width W of such a dent is given 
for t < T by 



1 + 


e o 


-i2 


V Z (T - t) 



(70) 


At great values for T - t, this expres- 
sion approximates 


W 3 [2V z e 0 (T - t )] 1/2 

= {2V z e 0 [(T - ti) - (t - ti )]} 1/2 

(7D 



Figure 9*- View of trace in |,z plane of 

Thus, with increasing t, W diminishes perturbed boundary with positive z,,!. 

with a rate dW/dt that increases with 

time. If, for example, at an initial instant t = ti, W = 1000 km, 

V z = 500 km/sec, e Q = 100 km, equation ( 71 ) indicates that T - tq = 10 sec. 

At a time 5 sec after ti, W = 707 km, and 5 sec later still equation (71) 
indicates that W = 0. Equation (71) is not sufficiently accurate for the 
latter time, however, and the more complete relation given by equation (70) 
shows, in fact, that W = e 0 . The initially shallow dent thus collapses to 
a semicircle in 10 sec. This configuration marks an obvious limit of which 
the above analysis may be employed to determine the shape of the entire dent. 
As indicated by equations (64) through (66), the center of the circle describ- 
ing the dent moves to positive z at subsequent times. If this solution were 
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used to describe the shape of the entire dent, part of the surface would be 
shielded from the solar wind, and the condition specified by equation (4) 
would not be satisfied. 

The appropriate solution for the model defined by equations (l) through 
(5) is as follows for the further development of the dent following the time 
t = T at which the dent is semicircular with the center of the origin. The 
portion of the dent that penetrates most deeply into the magnetosphere is 
semicircular and defined for z > z c by equation (64). The remainder of the 
boundary connecting the ends of the semicircle and the undeformed equilibrium 
portion of the boundary is straight and parallel to the z axis. This fol- 
lows from the fact that the magnetic pressure on this portion of the boundary, 
as indeed on the entire boundary in this example, is equal to B s 2 /&rt = KmnV z 2 . 
Since n s = ±i on these segments of the boundary, equation ( 3 ) shows that 

KmnV z 2 = Kmn(V x - v x ) 2 (72) 

Hence , 

|v x - Vx| = v z (73) 

where v x represents the velocity of the boundary, positive in the direction 
of the positive x axis, and the signs are associated with the two segments 
in such a way that the straight -sided intrusion into the magnetosphere col- 
lapses rather than expands with the passage of time. Viewed in the |,z 

coordinate system that drifts along 
the boundary with the tangential compo- 
nent of the solar wind V x , the 
straight segments travel toward the z 
axis with a uniform velocity equal to 
V z . The straight sides thus remain 
connected to the ends of the semi- 
circle, and all conditions are satis- 
fied. The resulting form of the 
indentation is illustrated in figure 10 
for the set of numerical conditions 
given above following equation ( 71 )* 

Of particular note is the result that 
the dent collapses from a semicircle 
to a line along the z axis in 
l/5 sec. 

Once a segment of the boundary 
becomes parallel to the z axis, it 
loses its ability to diminish the z 
component of the velocity of the inci- 
dent particles. Whether considered 
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from the classical Chapman -Ferraro point of view in which the individual par- 
ticles of the solar wind are specularly reflected from the boundary (K = 2 in 
eq. (3) ) , or from the more recently proposed point of view in which a fluid 
flows along the boundary (K = 1), particles encountering the straight segment 
of the boundary are merely shunted toward the opposite side of the dent as 
they continue their progress toward the circular nose of the dent. Although 
effects of a single particle experiencing multiple encounters with the bound- 
ary are not included in the foregoing analysis, it is evident that the number 
of particles, or the density of the fluid, confined within the dent increases 
with time as the dent collapses in width. As a result, the nose of the dent 
is subject to an additional pressure for a small interval of time immediately 
before the final collapse. This portion of the boundary will thereupon com- 
mence to move in the direction of the positive z axis, and it appears 
inevitable that a narrow column of compressed solar plasma will be injected 
in the magnetosphere. Its width will be a fraction of e Q , and quite unre- 
lated to the initial width of the dent. The time to collapse, however, 
depends very nearly on the square of the initial width, as indicated approxi- 
mately by equation (71). Further analysis and discussion of the behavior of 
the dent and possible effects of double or multiple encounters is provided in 
the following section. 


The actual magnetosphere boundary is not flat as assumed in this section, 
but curved. We have seen in the preceding section that curvature exerts a 
stabilizing or destabilizing influence on the boundary perturbations accord- 
ing to whether the center of curvature is on the field or plasma side of the 
boundary. For the former, equations (52) and (57) show that the amplitude of 
alined waves diminishes by a factor e in a time t C ur = R/V z . It is of 
interest to compare this time with the time collapse tcol of a shallow cir- 
cular dent of initial width W and amplitude e 0 * The ratio of these two 
times 


t R/V 2e R 

x cur ' z o 

tco1 W 2 /2V z e 0 w2 


( 7 ^) 


provides a measure of which tendency will tend to dominate. If t cur /t co q 
is much greater than unity, the stabilizing effect of curvature will be small, 
and a dent will tend to collapse much as in the case described above for the 
flat equilibrium surface. If, on the other hand, this ratio is much less 
than unity, the collapsing action will have little time to develop and the 
stabilizing effect of curvature will dominate. If, for a numerical example, 

R = 50,000 km and W = 1000 km, then t cur /t co q = 1 for e c = 10 km. Thus, 
such dents will tend to smooth out with time if e<o is substantially less 
than 10 km, and collapse with subsequent particle penetration if e Q is sub- 
stantially greater than 10 km. It would thus appear that even quite modest 
irregularities in the solar wind could result in isolated columns of solar 
plasma penetrating the boundary of the magnetosphere in spite of the stabiliz- 
ing influence of curvature. The plasma in these columns would, moreover, 
retain its tangential momentum as it initially penetrates into the magneto- 
sphere. Subsequent interaction with the ambient magnetosphere plasma would 
lead to a transfer of momentum across the boundary that might be roughly 
equivalent to the effective viscous interaction postulated in the 
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magnetosphere convection model of Axford and Hines (ref. 22). If, on the 
other hand, the solar wind is perfectly steady and free of irregularities, no 
penetration would occur, except in the vicinity of the neutral points where 
the boundary curvature reverses. 


MOTION OF PARTICLES DURING FINAL STAGES OF COLLAPSE OF 
A CIRCULAR DENT IN THE MAGNETOSPHERE BOUNDARY 


A critical point in the development of a semicircular cylindrical dent 
in the magnetosphere boundary is reached when the dent assumes such a form 
that a rebounding particle has positive rather than negative V z . Since a 
particle with positive V 2 will inevitably hit the boundary of the dent a 
second time, and add an effect not included in the foregoing analysis, it is 
of interest to inquire whether the dent reaches the semicircular form before 
the first double bounce occurs. If it does, the preceding discussion of the 
final stages of the collapse and possible injection of particles would have 
to be modified substantially. It is the object of this section to demon- 
strate, within the classical Chapman -Ferraro formulation of the problem with 
specular reflection at the boundary, that such double bounces do not occur 
before the dent becomes a semicircle. Double bounces do begin to occur 
slightly later, however, and the section is concluded with some quantitative 


In a £,z coordinate system 
drifting along the boundary with the 
tangential component of the solar wind 
and with the angles X and X 0 and 
other quantities as illustrated in fig 
ure 11, equations ( 3 ) and (5) combine 
to yield the following expressions for 
the normal velocity v n of a segment 
P of a circular dent 


V n /v z = -1 + cos X (75) 

Particles of the undisturbed incident 
solar plasma approach P with veloc- 
ity components V_l 0 and V- ^ perpendic 

ular and parallel the local surface 
element given by 

V_l q = V z cos X ) V|| = V z sin X 

(? 6 ) 


remarks on the consequences thereof. 
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Relative to P, which is itself moving with normal velocity v n , the corre- 
sponding component velocities of the particle are 


V_Lp = V z cos X - v n , V|| - V 7 sin X (77) 

n O 1 ! R 

n 0 

Upon impact, the tangential velocity component remains unchanged, and the nor- 
mal velocity component relative to P is reversed. Designating the velocity 
components after rebound by subscript 1, we have 


V_l r = -V z cos X + v n , Vn = V z sin X 
1 r i 


(78) 


Returning to the original |,z coordinate system, we have 


V x = -V z cos X + 2v n = V z (cos X - 2) , V,, = V z sin X 

1 I! i 


( 7 9) 


Resolution into components parallel to the |,z coordinate axes, that drift 
with the tangential component V x of the solar wind, yields 


V| = -Vj_ sin X + Vj | cos X = 2V Z sin X 


(80a) 


V z = Vj_ cos X + Vi| sin X = V z (l - 2 cos X) 


(80b) 


The resultant speed of the rehounding particles is thus (in the drifting coor- 
dinate system) 


1 / 2 


v x = (vf i + vf 1 ) =. V z (5 - h cos X 


1/2 


(8i) 


Equation (80h) shows that a particle that strikes the edge of the dent 
where X = Xo rehounds exactly along the | axis with V z = 0 if 

cos Xo = l/2. The critical value for Xo is thus 60°. Equation (8l) shows 
that the rehounding particle travels with velocity \/3V z and arrives at the 
origin at exactly the _s_ame instant as the center of the dent, since the half 
width of the dent is \l 3 times the distance from the origin to the center of 
the circular arc defining the dent when X D = 60° and the center is itself 
moving with velocity V z . It is thus evident that this particle will not 
experience a double collision prior to the time at which the dent attains a 
semicircular form. 
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Since it is not immediately evident that the particle to which attention 
is directed in the preceding paragraph is the first to actually experience a 
double collision, it is important to inquire further and determine the coor- 
dinates, at the moment the dent becomes semicircular, of all the particles 
that hit the dent in all preceding time. To be precise, consider the location 
at t = T of the particles that collided with any- part of the dent at an 
earlier time t c when the center of the dent is at 

z c = -V z At , = 0 (82) 

the radius R is 

R = e 0 + V z At ( 83 ) 

and the angle X G is 

X 0 = cos _ 1 (V z At/R) (84) 

where At = T - t c . The coordinates of an arbitrary point P at time t c 
are 


= -R sin X ( 85 a) 

r 0 

z p = R(cos X - cos X Q ) ( 85 b) 

where , as always, |x| < |Xo| • Particles colliding with the dent at t c thus 
travel for an interval of time At with velocity components given by equa- 
tion (80) and arrive at the following coordinates at T. 

| = |p c + At = R sin X(-l + 2 cos X 0 ) ( 86 a) 

z = Zp + V z At = R cos X(l - 2 cos X Q ) ( 86 b) 

r c 1 

This result shows that these particles are situated on the portion of a cir- 
cular arc of radius 


(I 2 + z 2 ) 1/2 = R 


1-2 cos X r 


= 


V z At 


(87) 


that extends over the range of angles given by 


x| < |x 0 


1 — 1 

H 

1 

c 

( 

c 


-1 V z At 

1 + e 0 /V z At 


C ° S e 0 + v z At 


( 88 ) 
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The results are plotted in dinens ionless 
form in figure 12. The various circular 
arcs define the location at t = T of 
the particles that collided with the 
dent at such earlier times that the 
ratio V z At/e 0 had the numerical value 
indicated. The final conclusion is that 
a circular dent in the magnetosphere 
boundary is able to collapse all the way 
to a semicircle before the first double 
bounce occurs. 

For a certain interval of time after 
the center of the circle passes the ori- 
gin, the dent will continue to collapse 
without any double encounter with a par- 
ticle. During this time, the nose of 
the dent will continue to be described 
by a semicircle as discussed previously 
and as illustrated in figure 12. The 
radius of the circular nose is 
e o " V z (t- T), and the sides are par- 
allel to the z axis. The flat sides 
of the dent are thus traveling toward 
the z axis with velocity equal to V z , 
and particles encountering these surfaces 
will be deflected with a velocity compo- 
nent in the | direction equal to ±2V Z 
while the component in the z direction 
remains unchanged and equal to V z . 

As demonstrated near the beginning 
of this section, the first particle that 
is deflected in such a way as to' insure 
a double impact is the one that makes 
initial contact with the edge of the 
dent when X Q = 60°. It is deflected 
so as to travel exactly along the | 
axis with velocity \/3Vz, and hence 
passes through the origin at the same 
instant (t = T) as the center of the 
dent. This particle continues on its 
way, and makes its second encounter with 
the surface of the dent at the instant 
when the radius of the circle is equal 
to e 0 (3 - \f 3)/2 » 0.63e o . The center 
of the circle is thus at about O.37eo> 
and the configuration of the dent is 
as illustrated in figure 13 . The addi- 
tional pressure exerted on the boundary 
by these and succeeding particles will 



Figure 12. - Location at time T of parti- 
cles that were specularly reflected from 
the dent at earlier time. 




Figure 13. - Configuration of dent at instant 
when particle which was deflected from 
edge of the dent when X 0 = 60° makes its 
second impact. 
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cause the portion of the flat surface nearest the axis to he slowed somewhat 
in its advance toward the z axis. This effect must he of a rather transient 
nature, however, because the initial deflection must occur after the time when 
X = 60°, and the entire effect is terminated well before the time when 



Figure 14.- Gotif Lguration of dent at instant 
when first particle to be deflected from 
the flat surface of the dent makes its 
second impact. 


X 0 = 90 . In any case, the total 
effect is expected to he too weak and 
transient to prevent this portion of 
the dent from continuing its advance 
toward the z axis until, eventually, 
the two sides meet and the original 
plane boundary of the geomagnetic field 
is restored along the £ axis. 

At a slightly later instant, when 
the center of the semicircle progresses 
to 2 e 0 / 3 , the first particle to bounce 
off the flat surface will encounter the 
opposite side of the dent. It will hit 
at exactly the junction of the straight 
and circular elements, as illustrated 
in figure 14, and moreover, will be 
followed by a steady stream of parti- 
cles extending over a steadily widen- 
ing portion of the side of the dent. 

The particles approaching their second 
encounter with the dent in x > 0 have 
velocity components V z = V z , = 2V Z , 
and number density n = n /2 where n 
is the density in the undisturbed 
stream. Then equations (3) and (5) 
combine to yield the following expres- 
sion for the velocity of the affected 
portion of the straight segment of the 
dent 


Bge/Stf = 2mnV^ = 2mn(V x - v n )‘ 


SEMICIRCULAR 



Z 


MULTIPLE- BOUNCE PENETRATION 
AND PINCHING OFF 





v 


Figure Sketch of possible mechanism for 

Injection of columns of solar wind, plasma 
in to magneto s phere . 


= mn(2V z - v n ) ( 89 ) 

Solving for v n /Vz> we get 

v n /V z = 2 - V2 ~ 0.59 (90) 

We thus find this portion of the dent 
bulging out into the magnetosphere at 
the same time as the portion of the 
dent near the £ axis is snapping 
shut much in the manner illustrated in 
figure 15. This analysis supports the 




more qualitative discussion provided in the preceding section, and lends sup- 
port to the conclusion presented there that this mechanism appears capable of 
injecting widely separated columns of solar plasma into the magnetosphere 
under all but the steadiest of conditions. 


Ames Research Center 

National Aeronautics and Space Administration 
Moffett Field, Calif., June 11, 1964 


REFERENCES 


1. Dungey, J. ¥. : Cosmic Electrodynamics. Cambridge University Press, 

Cambridge, 1958, pp. 151-152. 

2. Dungey, J. W. : The Structure of the Exosphere or Adventures in Velocity 

Space, Geophysics, the Earth's Environment. C. DeWitt, J. Hieblot, 
and A. Lebeau, eds., Gordon and Breach, Science Publishers, New York, 
London, 1963 , pp. 505-550. 

3. Parker, E. N. : Interaction of the Solar Wind With the Geomagnetic Field. 

Phys. Fluids, vol. 1, 1958, pp. I 7 I-I 87 . 

4. Hines, C. 0.: The Magnetopause: A New Frontier in Space. Science, 

vol. l4l, 1963, PP. 130-136. 

5. Spreiter, John R. : The Boundary of the Geomagnetic Field. ICSU Review 

of World Science, vol. 6, 1964, pp. 178-190. 

6. Chapman, Sydney: Solar Plasma, Geomagnetism and Aurora. Geophysics, the 

Earth's Environment. C. DeWitt, J. Hieblot, and A. Lebeau, eds., Gordon 
and Breach, Science Publishers, New York, London, 1963 , pp. 373-502. 

7. Dessler, A. J. : The Stability of the Interface Between the Solar Wind 

and the Geomagnetic Field. J. Geophys. Res., vol. 66, 1961 , PP- 3587- 
3590. 

8. Dessler, A. J. ; Further Comments on Stability of Interface Between Solar 

Wind and Geomagnetic Field. J. Geophys. Res., vol. 67 , 1962, pp. 4892- 
4894. 

9. Coleman, P. J., Jr., and Sonett, C. P. : Note on Hydromagnetic Propaga- 

tion and Geomagnetic Field Stability. J. Geophys. Res., vol. 66, 1961 , 
PP * 3591-3592. 

10. Barthel, J. R. , and Sowle, D. H. : A Mechanism of Injection of Solar 

Plasma Into the Magnetosphere. Planetary and Space Science, vol. 12, 
1964, pp. 209-217. 


29 


11 . 


Spreiter, John R. , and Alksne, Alberta Y. : On the Effect of a Ring Cur- 

rent on the Terminal Shape of the Geomagnetic Field. J. Geophys. Res., 
vol. 67, 1962, pp. 2193-2205. 

12. Spreiter, John R. , and Alksne, Alberta Y. : The Effect of a Ring Current 

on the Boundary of the Geomagnetic Field in. a Steady Solar Wind. 

NASA TR R-177, 1963. 

13. Spreiter, John R. , and Briggs, Benjamin R. : Theoretical Determination of 

the Form of the Hollow Produced in the Solar Corpuscular Stream by 
Interaction With the Magnetic Dipole Field of the Earth. NASA TR R-120, 
1961. 

14. Spreiter, John R. , and Briggs, Benjamin R. : Theoretical Determination of 

the Form of the Boundary of the Solar Corpuscular Stream Produced by 
Interaction With the Magnetic Dipole Field of the Earth. J. Geophys. 
Res., vol. 67, 1962, pp. 37-51. 

15. Axford, W. I. : The Interaction Between the Solar Wind and the Earth’s 

Magnetosphere. J. Geophys. Res., vol. 67, 19^2, pp. 3791-3796. 

16. Kellogg, P. J. : Flow of Plasma Around the Earth. J. Geophys. Res., 

vol. 67, 1962, pp. 3805-3811. 

17. Spreiter, John R. , and Jones, Wm. Prichard: On the Effect of a Weak 

Interplanetary Magnetic Field on the Interaction Between the Solar Wind 
and the Geomagnetic Field. J. Geophys. Res., vol. 68, 1963, pp. 3555- 
3564. 

18. Ferraro, V. C. A. : On the Theory of the First Phase of a Geomagnetic 

Storm: A New Illustrative Calculation Based on an Idealized (Plane not 
Cylindrical) Model Field Distribution. J. Geophys. Res., vol. 57, 

1952, pp. 15-49. 

19. Beard, David B. : The Interaction of the Terrestrial Magnetic Field With 

the Solar Corpuscular Radiation. J. Geophys. Res., vol. 65, i960, 

PP. 3559-3568. 

20. Mead, Gilbert D. : Deformation of the Geomagnetic Field by the Solar Wind. 

J. Geophys. Res., vol. 69, 1964, pp. 1181-1195. 

21. Fukushima, N. : Gross Character of Geomagnetic Disturbance During the 

International Geophysical Year and the Second Polar Year. Report of 
Ionosphere and Space Research in Japan, vol. 16, 1962, pp. 37-56. 

22. Axford, W. I., and Hines, C. 0.: A Unifying Theory of High-Latitude 

Geophysical Phenomena and Geomagnetic Storms. Canadian J. Phys., 
vol. 39, 1961, pp. 1433-1464. 


30 


NASA -Langley, 1965 A “995 



